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Abstract 

It has been argued that the speed-density digram of pedestrian movement has an inflection 
point. This inflection point was found empirically in investigations of closed-loop single-file 
pedestrian movement. The reduced complexity of single-file movement does not only allow a 
higher precision for the evaluation of empirical data, but it occasionally also allows analytical 
considerations for micosimulation models. In this way it will be shown that certain (common) 
variants of the Social Force Model (SFM) do not produce an inflection point in the speed-density 
diagram if infinitely many pedestrians contribute to the force computed for one pedestrian. We 
propose a modified Social Force Model that produces the inflection point. 

1 Introduction: Empirical Data on Pedestrians’ Speed- 
Density Relation 

In the course of recent years a number of experiments have been conducted in which pedestrians 
walk single-file in a closed loop [1-5]. Having a different number of pedestrians in the loop, 
different densities are prepared. In a section of the loop line density and speed are measured. 
Figure 1 shows the experimental setup from which most data stems. There were experiment 
runs in which the loop was larger and more pedestrians participated, but the principle was 
always the same. 



Figure 1: Experimental setup. Source: Figures 2 and 3 of [1] 

The experiment has been conducted at various places around the world. Figure 2 shows 
the results for India and Germany. With the free speeds plotted in this diagram the existence 
of an inflection point is obvious. It is not quite clear, however, at which density curvature is 
maximally negative and also the density of the inflection point can only be estimated roughly. 
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Figure 2: Speed density diagram for India (left) and Germany (right). Source: Figure 5 of [2] 


Whereas Figure 2 clearly shows that differences between India and Germany were found, the 
loop size (corridor length) and the profession of participants did not make a difference, at least 
not an obvious one. Compare Figure 7 of [2] with Figure 1 of [5]. 

The discussion 1 in [1] indicates that the existence of an inflection point is common sense and 
more so that for - at least moderately - high densities the curvature of the flow density relation 
is positive. 


2 The Social Force Model for Steady-States in Single-File 
Movement 


The circular specification of the Social Force Model [6] is defined as 2 : 
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where vg a is the desired speed of pedestrian a. A a > 0, B a > 0, 0 < X a < 1, and r a > 0 
are parameters of the model. R denotes the body radius of a pedestrian. e a p points from 
pedestrian (3 on pedestrian a. x is the position of a pedestrian and dots mark time derivatives 
(i.e. speed and acceleration). The sum runs over all - potentially infinitely many - pedestrians 
in a simulation scenario. Function w() is there to suppress forces acting from behind. In it 0 a p 
is the angle between pedestrian a's velocity vector and the line connecting a and /3. 

From here on we assume that parameters A, B , A, r I?, and vq have identical value for all 
pedestrians, so we omit the indices. This allows to combine A and R into a new parameter 
A = Ae~ 2R . 

Since it is obvious which properties are time dependent, we also omit the “(f)”. Then 
equation (1) reads for the one-dimensional case: 


^0 - x a y - 

x a = -1- A 2_^ uJ{x a , xp, X)e b ( 3 ) 

T P 

d a p \%p x a \ (4) 

w(x a ,xp, X) = Xiix/3—x a <0 (5) 

w(x a ,xp,X) = —lifxp—x a >0 ( 6 ) 


1 Quote (p. 3): “Domain I: ... At low densities there is a small and increasing decline of the velocity ... Domain 

III: ... For growing density the velocity remains nearly constant.” 


2 We neglect here forces from walls from the beginning since we are anyway working towards the one-dimensional 
case. 
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with the additional assumption that for all pedestrians and times x > 0. Thus we do not intend 
to model the loop from the experiment, but assume closed boundary conditions in Id. 

Elliptical specification II is a variant of the Social Force Model where the force between 
pedestrians - in addition to the mutual distance - also depends on the relative velocity of 
pedestrians a and /3. The full definition is give in [7] and reduced to one dimension it can be 
found in [8] . If - as in the scenario discussed in this work pedestrians a and /3 have the same 
velocity - vanishing relative velocity - elliptical specification II gives the same force for a as the 
circular specification. Therefore the further reasoning applies also to elliptical specification II 
(not elliptical specification I, though). 

Now we investigate the steady state of this in fact both - model(s). Steady-state means 
that the speeds and distances remain constant and consequently acceleration is zero for all 
pedestrians. With the left side of equation (3) being zero we can easily compute the steady 
state speed without having to solve a differential equation: 

x a = v 0 + tA^ w(x a ,Xff, X)e~~^ (7) 

P 


As we chose that the parameters for all pedestrians are equal and the system has periodic 
boundary conditions or is infinitely large, all distances between neighboring pedestrians are 
equal. Then necessarily the distance (center point to center point) from pedestrian a to the n 
next neighbor /3 n can then be written as: 


d a pn = nd 0 (8) 

If we resolve the w() function into both directions - w() = —1 for all pedestrians ahead and 
w() = A for all pedestrians behind - we can rewrite equation (7) more explicitly with the sum 
running over natural numbers not pedestrians: 


x a = v 0 - (1 - A)tA^ e B ° (9) 

n =1 

Since do and B both necessarily are positive it is e ^ <1 and the sum is the geometric 
series with the well known result 
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where p is the line density of pedestrians p = 1/do- 

With appropriately chosen values for parameters Vg. A, r A, B equation (12) should be able 
to reproduce the empirical fundamental diagram as shown in Figure 2 - obviously not each 
single data point, but the general, average course. This includes that the speed density relation 
in equation (12) computed for/from the Social Force Model should also yield an inflection point. 
So for (12) we compute the second derivative of the speed function with regard to density: 
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and test if it is zero for some density pc 


(2 Bp, - l)e^ - (2 B Pi + 1) = 0 


(16) 


This does not have a solution, but the left side of the equation is always negative and approaches 
zero only asymptotically for pi —>• oo. 

At this point we could write down conclusions and title the paper “Requiem for the Social 
Force Model” (and btw. - a number of other models as well). However, experience teaches 
to first search for possibilities of resurrection. Still we note as first remarkable result of this 
work: Result 1: Neither the circular specification nor the elliptical specification II of the Social 
Force Model as originally defined produce an inflection point in the speed density relation for 
homogeneous steady-state one-dimensional movement. 

The next step does not follow as a consequence of what is written here so far, but it has 
to be justified a posteriori: we reconsider equation (7) and investigate a variant of the model 
where not all pedestrians, but only the nearest neighbors - the one in front and the one at the 
rear - exert a force on pedestrian a. This leads to 

x a = vo — (1 — A )rAe~ Wp (17) 

which we write seemingly unnecessarily complicated as 

x a = v 0 - (1 - A)tA—- (18) 

— 0 


In this form we note that the only difference between equation (12) and equation(18) - so 
between the original formulation where infinitely many pedestrians exert a force on a and the 
nearest neighbor variant - is only that in the first case there is a “1” and in the second case a 
“0” in the denominator. 

Having written the equations in this way the natural next question is “What if instead of 
the zero or one in equations (12) and (18) we write there some 0 < k < 1?”: 
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gBp - 
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( 21 ) 
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(23) 


For 0 < k < 1 the new equation to solve is 

(2 Bpi — l)e Bp i — k{2Bpi + 1 ) = 0 


(24) 


and it has a solution. For k = 0 it is obviously pi = 1/(2 B). For other values of k the solution 
needs to be computed numerically. Table 1 gives some values for Bpi. 

The introduction of parameter k produces the desired inflection point in the fundamental 
diagram. See figure 3. To find an intuitively comprehensible interpretation for this extension, 
we write equation (21) in a slightly different manner 

v k (p) = v 0 - (1 - X)tAj^^ - (25) 

rZ e Bp 

~k 1 
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Table 1: Numerical solutions for the value pt of the inflection point with regard to various values 
for parameter k. Values computed with [9] 


k 

Bpi 

k 

Bpi 

k 

Bpi 

0.0 

0.500 

0.90 

0.981 

0.99 

2.049 

0.1 

0.515 

0.91 

1.013 

0.999 

4.379 

0.2 

0.531 

0.92 

1.051 

0.9999 

9.416 

0.3 

0.551 

0.93 

1.096 

0.99999 

20.28 

0.4 

0.576 

0.94 

1.151 

0.999999 

43.68 

0.5 

0.606 

0.95 

1.219 

0.9999999 

94.10 

0.6 

0.646 

0.96 

1.309 

0.99999999 

202.7 

0.7 

0.703 

0.97 

1.435 

0.999999999 

436.8 

0.8 

0.793 

0.98 

1.635 

0.9999999999 

941.0 


and undo the summation of the geometric series with k “on the back” of the exponential function 

, OO 

Vk(p) = Vq — (1 - A)ryl— k n e^T^ 

n—1 
oo 

= v 0 — (1 — A )tA ^ fc ri ~ 1 e~ g p 

n—1 

In terms of forces/acceleration this means: 

x = ^ - (1 - A )Ajr t k n - 1 e- S & (28) 

n—1 

This is easy to interpret: the next neighbors of a exert a force on a unmodified compared to the 
original SFM. The second next neighbors exert a force which is suppressed by a factor k (always 
compared to the original model without parameter k). The force from the next to next to next 
nearest neighbors is suppressed by a factor of k 2 and so on. This implies for example that if 
we take out each second pedestrian from a simulation and the next to next nearest neighbor 
becomes the next nearest neighbor the force from this new nearest neighbor is larger than before 
when s/he was just the next to next nearest neighbor although the distance to a is the same as 
before. 

This can be comprehended intuitively. If two pedestrians approach “me” independently and 
the first remote pedestrian overtakes and becomes the closest one “my” awareness is shifted 
from the former to the new nearest neighbor and that would have an impact on my velocity 
changes, i.e. acceleration, i.e. forces. This intuitive comprehensibility of the model extension 
adds to the pleasure of having gained the desired inflection point. Thus: 

Result 2: Suppressing the force from each pedestrian in the sequence ordered by distance from 
the pedestrian for whom forces are calculated with an additional factor 0 < k < 1 produces an 
inflection point in the resulting macroscopic speed-density relation and as a model extension can 
be motivated intuitively. 

We would like to emphasize that the introduction of parameter k as extension of the So¬ 
cial Force Model brings a major conceptual change. Without it forces of various pedestrians 
superpose without interfering. With parameter k on the contrary one has to know the local dis¬ 
tribution of all pedestrians before one can compute the force of one pedestrian on another one. 
Forces do not superpose anymore. Instead the extended model - let us briefly call it “SFMk” - 
is rather described by Sherif’s famous description of social systems: “the properties of any part 
are determined by its membership in the total functional system.” [10]. In this sense the SFMk 


(26) 

(27) 
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structurally bears some similarities to another extension of the Social Force Model - namely the 
dynamic potential [11-14] - where the desired direction of pedestrians is computed such that 
they walk into the direction of earliest expected arrival, for this considering the distribution and 
movement state of all pedestrians in a holistic way. A difference, however, between SFMk and 
the dynamic potential is that in SFMk it is exactly known which pedestrian /3j contributes what 
to the effect on the movement of pedestrian a. This is not the case for the dynamic potential. 
We would therefore call the dynamic potential a mesoscopic or mean-field or holistic modeling 
element while the SFMk is non-superposing but non-holistic and entirely microscopic. 




Figure 3: Speed-density (left) and flow-density (right) relations of the /c-extended Social Force 
Model for various values of parameter k. The plots are based on the dimensionless variants of 
equations and it is a = 0.2. 


3 Summary 

In this contribution we have computed the macroscopic limit of the Social Force Model of 
pedestrian dynamics for single-file uni-directional movement. We found that the speed-density 
curve does not have an inflection point, but that with the introduction of one more parameter 
to the Social Force Model one can recover an inflection point. The new model intuitively 
makes sense considering that humans focus their awareness. 
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